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ABSTRACT 

The splitting of photons into two photons becomes both possible and significant in 
magnetic fields in excess of 10 12 Gauss. Below the threshold energy 2 m r c 2 for single- 
photon pair production, splitting can be an astrophysically-observable phenomenon evident 
in gamma-ray burst spectra. In such circumstances, it has been found that magnetic pho- 
ton splitting reprocesses the gamma-ray burst continuum by degrading the photon energy, 
with a net effect that is quite similar to pair cascade reprocessing of the spectrum. Results 
are presented for the spectral modifications due to splitting, taking into account the differ- 
ent probabilities for splitting for different polarization modes. Unpolarized and polarized 
pair cascade photon spectra form the input spectra for the model, which calculates the 
resulting splitting-reprocessed spectra numerically by solving the photon kinetic equations 
for each polarization mode. This inclusion of photon polarizations is found to not alter 
previous predictions that splitting produce a significant flattening of the hard X-ray con- 
tinuum and a bump at MeV energies below a pair production turnover. The spectrum 
near the bump is always strongly polarized. The spectra of bursts detected by GRO may 
provide observational confirmation of magnetic photon splitting. 

INTRODUCTION 

Magnetic photon splitting 7 — > 77 is an exotic and comparatively recent prediction 
of quantum electrodynamics in strong magnetic fields: the first correct calculations of the 
reaction rate were performed by Bialynicka-Birula and Bialynicki-Birula (1970), Adler et ai 
(1970) and Adler (1971). When B ~ 0 , splitting is kinematically possible but is forbidden 
by a charge conjugation symmetry of QED known as Furry’s theorem (e.g. see Jauch 
and Rohrlich, 1980). The theorem states that ring diagrams that have only external lines 
corresponding to photons, and that have an odd number of vertices, contribute zero to the 
interaction matrix elements. This symmetry is broken by the presence of an external field. 
The splitting of photons is therefore a purely quantum effect, and has appreciable reaction 
rates only when the magnetic field is at least a significant fraction of the quantum critical 
field B c — m^c 3 /(eh) — 4.413 x 10 13 Gauss. Neutron stars, and especially gamma-ray 
bursts (GRB) associated with neutron stars, with their high 7 -ray luminosities, therefore 
provide good environments for testing for the existence of splitting. 

The potential importance of magnetic photon splitting, suggested by Baring (1988) 
and Mitrofanov et al. (1986), has been explored in detail by Baring (1991) in the context 
of spectral implications for gamma-ray bursts. Splitting was found to degrade the photon 
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continuum in a manner similar to the reprocessing effect of a pair cascade. In magnetized 
neutron star environments, where B 10 12 Gauss, 7 — > 77 produces a flattening of 
the hard X-ray continuum and a broad bump at energies around 1 MeV (Baring, 1991). 
This feature is potentially observable, and if detected, would be a significant verification 
of QED in strong magnetic fields. In particular, the instruments aboard the GRO mission 
have the capability to make such a detection. In the model of Baring (1991), the source 
of continuum photons was generated by a synchrotron radiation pair cascade initiated by 
some quasi-isotropic injection of non-thermal electrons. 

Here the basic nature of the splitting process is reviewed, and the spectral modi- 
fications to the GRB continuum reconsidered: the new development in this work is the 
inclusion of an exact treatment of the individual polarization states of the photons par- 
ticipating in the splitting process. This extension of the work in Baring (1991) is an 
important step, particularly since it has long been recognized that splitting may be an 
important polarizing mechanism in neutron star sources such as pulsars (Adler et a/., 
1970; Usov and Shabad, 1983). It is found here that consideration of the photon polariza- 
tion states does not qualitatively change the spectral reprocessing nature of 7 77 , and 

that the MeV bump still appears. However, the resultant spectra are strongly polarized 
by the splitting process, and any detailed GRB spectral model should include the photon 
polarizations in order to make quantitative predictions even for observational detections 
without polarimetric capability. 

In this paper, it is convenient to use a simple dimensionless notation for various 
quantities in relativistic, magnetized environments. Therefore all magnetic fields will be 
expressed hereafter in units of the critical field B c , and all photon and particle energies 
will be scaled by m e c 2 . 

MAGNETIC PHOTON SPLITTING 

Photon splitting is a third-order QED process that must have a reaction rate that is 
strongly dependent on the magnetic field strength. The calculation of its rate by standard 
QED techniques is both difficult and laborious, because of the complications introduced 
by the presence of the field. The most formidable difficulty is the complicated form for 
the electron propagator in the external field, which renders the evaluation of the general 
expression for the splitting rate analytically intractable. The earliest research on this 
process therefore specialized to the low photon energy e<l (in units of m r c 2 ), low field 
(flCl) limit where dramatic simplifications can be made. This limit corresponds to 
the case of small dispersion. Since the presence of a magnetic field polarizes the vacuum, 
dispersive effects naturally occur. However when the field is low ( B <> 0.1 ) and the 
energy of photons that are splitting is moderate ( t £ 1 ) these dispersive effects are small: 
this corresponds to the situation in gamma-ray bursts (Baring, 1991). The most general 
calculation of splitting valid for any values of e and B was performed by Stoneham (1979). 

The calculations of Bialynicka-Birula and Bialynicki-Birula (1970), Adler et al. 
(1970) and Adler (1971) were all in this limit of low dispersion. An expression for the 



splitting rate for this case, averaged over all photon polarizations, was obtained by Papa- 
nyan and Ritus (1972). For the splitting of photons of energy em v c 2 with e 1 into 
photons of energies u>m r c 2 and (e — u>)m e c 2 , their result can be expressed as an optical 
depth for an emission region of size R (see Baring, 1991): 
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Here af ^ 1/137 is the fine structure constant, A c = h/{m e c) is the Compton wavelength, 
and 0 is the angle of propagation of the photon that splits with respect to the field 
lines. This is also the propagation angle for both of the produced photons: by energy 
and momentum conservation the splitting process is collinear in the absence of dispersive 
effects. The total rate for photon splitting is obtained by integrating eq. (1) over , 
yielding 



For neutron star atmospheres, r 8p (e) can be greater than unity if 0 is not too small. 

Adler (1971) observed that in the low-energy limit, the splitting rate had the same 
field dependence but a different magnitude for different combinations of the initial and 
final photon polarizations. Therefore, photon splitting must be a powerful polarizing 
mechanism. The polarization-dependent rates can be taken from eq. (23) of Adler (1971), 
which can be related to eq. (1), using the optical depth notation preferred here, via 
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The photon polarization labelling convention of Stoneham (1979) is adopted here (this 
standard form was not used by Adler, 1971): the label || refers to the state with the 
photon’s electric field vector parallel to the plane containing the magnetic field and the 
photon’s momentum vector, while i_ denotes the photon’s electric field vector being normal 

to this plane. The splitting modes _L— >_L|| , || — > J L and ||— >|||| are forbidden in the low 

energy, low dispersion limit. 

In order for photon splitting to have observable effects in gamma-ray burst spectra, 
not only must its optical depth in the GRB emission region exceed unity, but it must not 
also be dominated by any other photon absorption process. Its main rival as a mechanism 
for gamma-ray attenuation is single-photon pair production 7 — > e + e~ , which has a 
kinematic threshold at e sin 9 = 2 . Above this threshold, which is strongly dependent on 
9 , the reaction rate for 7 — > e^e~ is so great that it totally dominates photon splitting and 
renders the GRB emission region opaque to gamma rays (e.g. Baring, 1989). Therefore, as 
was deduced by Baring (1991), there is a small but significant energy range below the pair 
production threshold for which photon splitting can act effectively to modify the gamma- 
ray spectra, provided that 0 is not small. This is sufficient for photon splitting to create 
observable effects in the MeV continuum of burst sources. 
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SPLITTING EFFECTS IN GAMMA-RAY BURST SPECTRA 

The principal effect of photon splitting on a gamma-ray spectrum is to reprocess the 
spectrum by degrading the energies of the photons. A pile-up of photons might therefore 
be expected at energies near where the splitting optical depth in eq. (2) drops below unity, 
which is about 0.5 MeV when 0 ~ 7r j 2 and B ~ 0.1 . To determine the modifications that 
photon splitting makes to the gamma-ray burst continuum, it is appropriate to numeri- 
cally solve the time-independent kinetic equations for the two photon polarization states. 
The solution then has two equilibrium populations, nj_(e) and TZ|j (e ) , for the different 
polarizations. In the simplest model, photons can be injected into the neutron star mag- 
netosphere with some energy distribution g(e) for a specific emission angle 0 and then be 
allowed to split or escape from the emission region (the effect of pair production can be 
included in q(e) ). Since photon splitting is a collinear process, the kinetic equation can 
be solved separately for each 0 , and then integrated over emission angles. 

The kinetic equations can be written down quickly: the injection rates (possibly 

polarized) are q± ) \\{€)/t vsc and the photon escape rates are n^ \\(e)/t VHC , where fos c R/c 

is the photon escape time for the emission region. Note that here spatial and angular 
diffusion of the photons due to splitting is negligible because 7 — » 77 is an angle-preserving 
process. Since photon splitting absorbs photons and also creates them at lower energies, it 
provides both source and loss terms to the kinetic equations. If cr denotes the polarization 
modes _L and || , the loss terms are the absorption rates r fT (e)n rT (e)/t rsc ( T a (e) is defined 
in eq. ( 4 b) ), while the source terms are proportional to the integral of the spectral rates 
in eq. (3) when weighted by the photon populations n (r (e) . The kinetic equations are 
therefore integral equations. Scaling by X sr , these equations assume the form 
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is the absorption optical depth for polarizations cr =X , || , obtained using eq. (3). Note 
that the term with the integral includes an extra factor of two to account for the duplicity 
of photons created in the splitting process. If the photon polarizations are unobservable, 
these two integro-differential equations reduce to one. 

The injection spectrum q(e) is chosen to be a typical output of the synchrotron 
pair cascade model of Baring (1989). The continuum spectrum is therefore formed by 
synchrotron radiation, though it could also be produced by other mechanisms such as 
inverse Compton scattering and curvature radiation. For synchrotron photons supplied by 
non-thermal electrons, generally a power-law continuum results, though this is truncated 
due to the strong absorption of gamma-rays by pair production 7 — > e + e - (Baring, 1989). 
A canonical spectrum representing the injection expected from a synchrotron pair cascade 
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is then 
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where the optical depth of the emission region (of size R ) to single-photon pair production 
can be obtained from Baring (1991): 
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This approximation is valid for 2 < e± 1/5 , i.e. near the threshold for pair production. 
The spectral index a of the injection in eq. (5) is independent of the polarization state since 
synchrotron radiation and most other radiation mechanisms do not produce polarization- 
dependent indices from non-thermal electrons. Validity of an injection such as in eq. (5) 
requires that the pair cascade spectrum is generated much faster than the photons can 
split. This injection is therefore physically reasonable because synchrotron radiation and 
pair production generally have much shorter timescales than does photon splitting. 

The solutions of eq. (4) for an unpolarized injection of photons from a pair cascade 
are displayed in Fig. 1 . The reprocessing effect of magnetic photon splitting generates an 
absorption of photons just below the pair production threshold and an overproduction of 
photons at lower energies. The net effect is a flattening of the hard X-ray continuum and 
a broad spectral bump around 0.5 MeV when the photon emission angle is 0 — 7 r /2 . This 
bump appears in both polarization states and is qualitatively similar to the prediction 
made in the work of Baring (1991), where consideration of the photon polarizations was 
neglected. It is clearly an observable feature, and is moved to higher energies and becomes 
less prominent as 6 decreases. The pair production truncations to the hard gamma-ray 
spectrum are evident in Fig. 1 and generally occur at energies e ~ 2 / sin 6 for B ~ 0.1 . 
The polarization V(e) of the photons at energy e is given by 
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and is clearly quite large in Fig. 1 . In fact, in the idealized situation that the emission 
region has an enormous optical depth to j — > 77 , repeated splitting of photons generates 
a strongly-polarized photon population, 5/7 of which is in the JL state: this result can be 
obtained from a close consideration of eqs. (3) and (4). This preferential polarization of 
the continuum, corresponding to V = 3/7, is realized near the 7 — > e+e - cutoff for the 
sin 9 = 1.0 case in Fig. 1 . 

The photon injection from a synchrotron pair cascade will not be unpolarized, as 
assumed in Fig. 1 , so a more realistic injection is needed. Standard synchrotron radiation 
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Fig. 1: The solutions of the photon kinetic equations in eq. (4), illustrating the repro- 

cessing effects of splitting for different emission angles 6 . The polarization components 
n _i_( e ) (solid curves) and n||(e) (dashed curves) of the photon spectrum are shown. The 
injected photon spectrum is a truncated power-law as in eq. (5) with a = 0 , and is 
unpolarized (Q_l = Q|| )• The field is B = 0.1 . The dotted curves depict the polariza- 
tion V = |(nx — Ti||)/(n_L + n.||)|. Reducing 9 clearly diminishes the reprocessing and 
polarization effects of splitting. 

theory for power-law electrons predicts that a synchrotron spectrum of index a will have 
a corresponding polarization V = (a + l)/(a + 5/3) (e.g. see Bekefi, 1966), independent 
of the photon energy. This corresponds to nx(e)/ n ||( e ) = 3a + 4 , which for a = 0 is 4 . 
The pair production process only generates non-thermal pairs and so does not influence 
the synchrotron polarization. Fig. 2 shows the photon splitting-reprocessed spectrum fot 
the case of a synchrotron injection spectrum that is polarized according to the above pre- 
scription (i.e. Q i_jQ\\ = 3 a + 4 ). Clearly the synchrotron polarization tends to dominate 
the emergent photon polarization, without affecting the prominence of the splitting bump. 
However, in contrast to Fig. 1, near the pair production cutoff, where photon splitting is 
most effective, 7 — » 77 actually acts to depolarize the synchrotron radiation towards the 
value of V = 3/7 that is “preferred” by magnetic splitting. 

Modelling gamma-ray burst spectra could plausibly include an integration over e- 
mission angles 6 : indeed observations are unlikely to sample just a single value of 6 . 
This can be represented by an enveloping of spectra like those in Fig. 1, resulting in Fig. 
3 when summed over polarizations. The angular distribution of radiation is chosen for 
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Fig. 2: The solutions of the photon kinetic equations in eq. (4), for a case almost identical 

to that depicted in Fig. 1: here the photon injection is again a truncated power-law with 
a = 0, but this time it is polarized as would be expected from synchrotron radiation (i.e. 
Q±/Q || = 3a + 4 — 4). The reprocessing due to splitting is again evident, however now 
the polarization is mostly determined by the synchrotron process. 

convenience (see Baring, 1991) to be (f)(s\n0) = sin* 6 , and negative values of k occur in 
models where the emitting electrons cool in energy and decrease their pitch angle during 
their synchrotron emission lifetime (this amounts to a pile-up of electrons with small pitch 
angles). The high-energy spectral break is by an index (see Baring, 1991) of (1 + a + k) , 
which becomes zero when k = —(1 + a) . It is due to the enveloping of the pair production 
cutoffs for different 9 . The photon splitting bump in the continuum around 1 MeV, which 
is still polarized, remains very prominent. The bump appears near 1 MeV because the 
dominant contribution to the splitting reprocessing comes from 6 ^ 7 t/ 2: it is unlikely to 
be confused with a two-photon pair annihilation line because it is so broad (Baring, 1991). 

In conclusion, magnetic photon splitting can produce a significant bump in the 
gamma-ray spectrum of GRBs, a prediction that is not altered by the accurate treat- 
ment of the photon polarization states that is presented here. The instruments aboard the 
GRO mission have the capability to detect this feature in burst sources, which would then 
be an important verification of the theory of quantum electrodynamics in strong magnetic 
fields. 
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Fig. 3 : Photon spectra resulting from an integration over emission angles and a sum 

over photon polarizations of spectra like those in Fig. 1, for B = 0.1 . The sum is weighted 
by the angular distribution <£(sin0) oc sin K 8 , with k as labelled. The photon injections 
are unpolarized ( — Qj| ), truncated power-laws as in eq. (5), with a — 0 (solid lines) 

and a = 1 (dashed lines). The splitting bump is more pronounced for flatter spectra, and 
the overall spectral break is (1 + a + k) . 
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